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Abstract 


The intrinsic high frequency permeability spectra of ferromagnetic composites 
containing different volume fractions of iron and cobalt have been simulated. A law 
(called Mghan’s law) is proposed to explain the simulated results by assuming that there 
are plenty of LLG (Landau-Lifshitz-Gilbert) type natural resonances contributing to the 
intrinsic permeability spectra. The results clearly show that the spectra strongly depend 
on the distribution of local effective magnetic field, the interaction between the magnetic 
particles, the inhomogeneous damping constant of LLG precession, and the initial 
equilibrium states. Especially, the effect of particles shape distribution in each sampling 
on the local effective magnetic field. In view of this fact: it is absolutely impossible to 
have the same effect from these factors when someone prepares several measurement 
samples, an uncertainty principle is believed to hold for measuring the intrinsic 
permeability of an electromagnetic (EM) composite. Therefore, this law tells us that it 
should be cautious when comparing or evaluating the EM properties of composites (for 
instance, EM wave absorbing composites). Memory effect can be used to restore the 


intrinsic high frequency permeability for a specific defunct composite sample. 
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1. Introduction 


Understanding the magnetic permeability (w) at high frequency is critical for developing 
magnetic devices or magnetic materials working at high frequencies. Electromagnetic attenuation 
requires larger magnetic imaginary parts (w") of complex permeability for dissipating the 
electromagnetic energy via large magnetic losses [1]. The planar magnetic inductors or magnetic 
microwave devices demandy” to be smaller for lower power losses [2]. It was also proposed that the 
imaginary parts of permeability can be negative (w"< 0) under the effect of spin transfer torque [3]. 
There are two well-known laws governing the relationships between permeability and ferromagnetic 
resonance frequency. One is “Snoek’s law” which describes the inverse relationship between the 
initial permeability and the loss peak of permeability spectrum [4]. The other one is the Acher’s law 
[5]: a revised Snoek’s law in an integral form of the product of (u” xf), which states that the integration 
of a magnetic loss spectrum must be smaller or equal to a constant (4nMs, Ms: saturation 
magnetization). However, both laws do not tell us the relationship between the frequency dependence 
of intrinsic permeability and microstructures (also includes magnetic domain structures) of magnetic 
materials, and the underlying physics governing the intrinsic permeability of magnetic composites. 
Especially, for an electromagnetic composite, it contains many ferromagnetic particles (or ferrite 
particles) which are in the unsaturated states, with many irregular shapes and are randomly distributed. 
These two laws are not applicable for explaining the variation of intrinsic permeability dispersion. 
Furthermore, in many published journal papers, a lot of electromagnet composites were claimed to 
have excellent electromagnetic properties, for instance, the electromagnetic wave absorbing 
properties. Can these “excellent” properties be reproduced experimentally by others? Or in other 
words, how much chances are there to reproduce their results? Both the Snoek’s law and Acher’s law 
also cannot answer these questions. In this paper, we will delve into these issues. In addition, it is 
frequently observed in practices that the spectrum of permeability (4z’(f)) has a broad peak, especially 
for a magnetic composite with inhomogeneous magnetic structures and material microstructures [6,7]. 
For magnetic composites, there are many proposed effective medium laws to retrieve the intrinsic 
permeability spectra from the data of composites made of magnetic particles (flakes commonly seen) 
and non magnetic matrix (for instance, polymers) [6]. Usually, a broad intrinsic permeability 
spectrum can be obtained. Here, the “intrinsic” permeability means that the permeability wholly 
attributed to magnetic materials excluding the dilution influences from non magnetic matrix, and the 
eddy current effect. More questions on this intrinsic permeability are required to be answered. What 
factors determine the dispersion behaviors of permeability, in other words, why the permeability 
values significantly vary with frequency? What factors determine the shape of the spectrum? How 
the interactions between the magnetic particles strongly impact on the obtained the intrinsic 
permeability? We believe that the intrinsic permeability of magnetic components significantly 
impacts on the measured effective permeability of composites. Therefore, understanding the 


frequency dependence of intrinsic permeability will be critical and helpful for developing the high 
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frequency magnetic devices or composites. To answer these interesting questions, in this paper, we 
are going to employ a micromagnetics simulation tool widely used in magnetism society to study the 


intrinsic high frequency permeability of a system with two magnetic phases. 


2. Simulation backgrounds and details 

Micromagnetics is a continuum theory to study the static behaviors of magnetization (hysteresis, 
magnetic domains) or dynamic responses of magnetization. Micromagnetics describe magnetization 
on a significant length scale which is large enough to replace atomic magnetic moments by a 
continuous function of position and small enough to reveal the transitions between magnetic domains. 
The micromagnetics simulations are performed using a three-dimensional object-oriented 
micromagnetics framework (OQOMMEF) by solving the Landau-Lifshitz-Gilbert (LLG) equation as a 


function of time [8-9] 


™ = —y(M x Herp) + < (Mi x ZA (1) 


dt 
, where Ms is the saturation magnetization and Hey is the effective field taking into account the 
exchange, demagnetization field, anisotropy and applied field terms. yis the Gilbert gyromagnetic 
ratio (2.21 x10° m/A.S), and @ is the damping constant for magnetization precession. In this work, 
the frequency dependence of permeability of a ferromagnetic plate (400 nm x 200 nm x20 nm, i.e. 
Length(L)x Width(w)xThickness(T)) are simulated by following these procedures. Firstly, the 
equilibrium configuration of magnetization is acquired in the absence of external applied magnetic 
field. Secondly, a weak pulse magnetic field with a form of h(t) =100exp(-10°) (t in s, h in A/m) is 
applied along the “z” axis (see Fig. 2), and h(t) is weak and is unable to drive magnetic domain walls 


to move. The temporal response of magnetization can be expressed as: 


6M (ri, t) =M(7,t) — Mort) = xr tht)  ---(2) 
The dynamic response of magnetization over the whole volume of sample is recorded. Both the pulse 
field and excited magnetization are then processed by a fast Fourier transform (FFT) approach, after 


which the relative permeability spectrum is calculated as below. 
1 r " SF 1 # " 
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, where fis the frequency, m(f) and h(f) are the expressions in frequency domain for M and pulse 
field (h) after FFT treatment, respectively; y ‘refers to the real parts of susceptibility and y” denotes 
the imaginary parts of susceptibility. The default parameters in OOMMF for Co phase are used in 
our simulations: saturation magnetization (Ms ) is 14x10° A/m, exchange stiffness constant (A) is 
30x 10°! J/m, anisotropy constant (K7) is 5.2x10° J/m?. For Fe phase, M; is 17 x105 A/m. A is 21x10" 


12 Ki is 4.7x10* J/m?. For all the simulation. The simulation object is discretized into many 
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tetrahedrons with a cell size of 2.5 nm x 2.5 nm x 2.5 nm. The analytic solutions of u'(f) and w"(f) 
of natural resonance from the above LLG equation can be expressed as follows [10], in which “œ” is 
the angular frequency (@ = 2rf): 


2 
wl = 145 [My?Ho(y?Hg — w?) + 20?) — (6) 
0 


" 1 
u" =+ aw(y?H3 + w?) —(7) 


D = (y2H2 — w?)? + 40? S ---(8) 
In these equations, Ho can be expressed as follows, which approximately equals to Herr since a (0.5 
used in this paper ) is greatly smaller than (yM). f- is the resonance frequency of one LLG natural 
resonance. 
az 
y? M2 


2 
Ho = iy, = Her p(1 + )/2 ~ Hepp ---(9) 


3. Results and discussions 


3.1. intrinsic high frequency permeability of two magnetic phases 


In our simulations, a nanostructured magnetic system is set up with two magnetic phases with 
different volume fractions using the method called “Voronoi diagram” [11]. As shown in Fig. 1. The 
Fe and Co nanoparticles inside a ferromagnetic plate (simulation object) are randomly oriented and 
in irregular shapes, which is alike in the microstructures of actual magnetic composites, for instance, 
thin films prepared by a co-deposition technique. The volume fraction of iron phase is 90.37 % and 


62.49 % respectively in Fig. 1 a, b. 


Co 


Figure 1. Microstructures of composites (a) Fe 90.37% + Co 9.63%; (b) Fe 62.49% + Co 37.51% 
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The initial spontaneously magnetized states for the case with both phases coexisting and for the case 
with only one Fe phase (or only Co phase) are compared, as shown in Figs. 2,3 and 5. In Fig. 2, the 
magnetization distributions in each composite sample are inhomogeneous, which are clearly shown 
by the colors and directions of arrows. The equilibrium direction of a magnetization vector depends 
on the local effective magnetic field (Hey). It is indicated that the Heg values within the sample are 


different and have a wide distribution. 
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Figure 2. (a) Fe 90.37% + Co 9.63%; (b) Fe 62.49% + Co 37.51% 

It is believed that it will influence the frequency dependence of relative permeability (called 
“permeability dispersion” in this work). For some applications (e.g., planar magnetic inductors) [12], 
the magnetic materials were magnetized to saturation so as to make almost all the magnetization 
pointing to the same direction. In this case, the permeability dispersion spectrum showed a sharp 
shape with high p’ values (real parts of complex permeability). Especially, the low u” values (ac 
magnetic losses) were obtained which was one of technically essential requirements. The narrow 
magnetic losses spectra w"(f) with the loss peaks being away from the working frequency are 
necessary for many electromagnetic devices. On the other hand, some applications require large high 
frequency magnetic losses, such as electromagnetic noises suppressors [13,14]. In these applications, 
magnetic composites contain plenty of magnetic particles which are usually not magnetized to 
saturation, or just in their spontaneously magnetized states, which is similar to conditions shown in 
Fig. 2. The fact that the local magnetizations are randomly oriented means there is a wide distribution 
of local effective fields, and results in the high-frequency permeability spectra spreading out, which 


will be shown soon. 
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Fig. 3. Composites only considering the different volume fractions of Fe particles (a) Fe, 


vol.90.37%; (b) Fe, vol. 62.49%. 


In order to show the law governing the intrinsic high frequency permeability of magnetic materials, 
especially with two ferromagnetic phases and with many particles in irregular shapes, we decompose 
the two phases of the composite plate into two parts without changing their original positions of 
particles in the following separate simulations. These two parts resemble jigsaw puzzles, which are 
important to eliminate the effects of altered particle position and particle size distribution. It is critical 
that if this prerequisite is not abided by, then the composite will not be the same, and the discussions 
will be meaningless. The equilibrium M configuration of only Fe phase in composites with different 
volume fractions are illustrated in Fig. 3, where the empty zones are the cobalt particles are originally 
located (where Co phase: A = 0; Ms = 0; K = 0). Subsequently, the permeability spectra of each phase 
are simulated separately for composites with different volume fractions, as shown in Fig. 4 (only for 
Fe phase) and Fig. 6 (only for Co phase). In both cases, the permeability spectra have a feature of 
spreading out. For most of electromagnetic applications (including this work), the microwave 
magnetic field is so weak and unable to drive the domain wall (DW) to move. There is no contribution 
from DW resonance mechanism (usually with f-_pw in MHz) to the permeability spectra. Natural 
resonance mechanism is the only one, as shown in LLG equation (1). However, the LLG equation 
tells us the behavior of natural resonance of only one magnetization vector. LLG natural resonance 
also works for the magnetic materials in a magnetized saturation state (single domain): we can treat 
the saturated magnetization as one “macrospin” [12]. However, the electromagnetic composites 
containing a lot of magnetic fillers which are not in the saturated state. Generally, their permeability 


dispersion spectra are very similar to the ones shown in Fig. 4 [7]. 
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Figure 4. Spectra of composites only considering Fe phase with different volume fractions.(a) Fe, 


vol. 90.37%. (b) Fe, vol. 62.49%. 


Here, it is assumed that the observed broad spectra are resulted from the coexistence of many LLG 
natural resonances. For one of initial equilibrium magnetized states (1,2,3...), a generalized law 
(called Mghan’s law) is proposed that summing up the permeability dispersion spectrum of each LLG 
resonance will result in the observed spectra. Suppose the magnetic material is in one magnetic 


66 9, 


equilibrium state “s”: 


at "s" state: u(f) = È iwi Xu, a)i) + € ---(10) 


, where: 


Vlei 
fo =a 11) 


27 


= = = 


Hey = Hk + Hs ---(12) 


— Kso 
D= nu 


1 
Ksq) =Z Ho X AN; xX MS ---(14) 


, where w; is the weight of one LLG natural resonance (i) in the entire resonances (n). AN is the 


demagnetization factor difference between the hard axis and the easy axis. For each one LLG 
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resonance, it has a specific dispersion of permeability < p(f,a); >, which has a resonance frequency 
J(i) and a specific damping constant aj. Hei) is the local effective anisotropic field. As we know, 
there are a lot of irregular particles (such as flakes) in a composite, AN obeys the lognormal 
distribution therefore, so do the Hei) and fxi) values according the above equations. €is an error term 
due to the interactions between magnetic phases and nonuniform damping constants, which will be 
shown soon. From the perspective of engineering practices with initial spontaneously magnetized 
states, the distribution of Hx and Hs should be paid more attentions. Hx strongly depends on the 
selection of materials (i.e., compositions). Hs strongly depends on the fabrication techniques (such 
as ball milling, thin film deposition, microwires drawing). The distribution of Hs (thus u(f) spectra) 
can be controlled by generally using a sieve to control the lognormal distribution parameters (such 
as “mode”, “median”, “mean”, or “standard deviation”) of AN for a pile of magnetic powder. If 
attempting to control the permeability spectra by controlling the remanence states (it will be discussed 
soon), the distribution of Hex also plays an important role and should be included into Eq. (12). The 
mathematical expression for the distribution of Hex is very complicate, but it can be mapped by FFT 
technique, evidences were given in our published paper [15]. As per one LLG natural resonance, the 
resonance frequency (f-) of Fe phase is calculated to be 1.58 GHz if only the magnetocrystalline 
anisotropic field (Hp) is considered (Hx=2k/(uoMs)). Taking into account the aspect ratio (L/T) of the 
Fe particles, the shape anisotropic constant (Ks) can be expressed as: K; = 0.5x1oxANxM,2. AN is 
roughly to be 0.5 for Fe particles in the composites. The f- value is found to be 31.0 GHz if 
simultaneously considering the magnetocrystalline and shape anisotropic field (Hs): (He = Hr + Hs). 


f- (31.0 GHz) is around the peak frequencies simulated and shown in Fig. 4. 


Figure 5. Composites only considering the different volume fractions of Co particles (a) Co, vol. 
9.63%; (b) Co, vol. 37.51%. 
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Next, let us see the results if the composites only containing cobalt particles are simulated. The 
equilibrium states of magnetization of Co phase are illustrated in Fig. 5, where the empty areas (set 
Fe phase: A = 0; Ms = 0; K = 0) represent the Fe particles located in the original nanocomposites in 
Fig. 2. Similar inhomogeneous distribution of magnetization vectors indicates the inhomogeneous 
distribution of local effective magnetic field. The permeability spectra are also expected to spread 
out, which are simulated and shown in Fig. 6. We also can explain the shape of permeability spectra 
using the proposed Mghan’s law shown in Equations 10-14. If only Hx is considered, a typical LLG 
natural resonance is to have the f, value as 20.8 GHz. If both Hy and H; are taken into account, f, 
should be 45.4 GHz, which can not be found in Fig. 6a. According to Mghan’s law, it is believed that 
both H, and Hex have a wide distribution that result in a wide distribution of f, values [15]. When 
there are no Fe particles interacting on the Co particles in Fig. 5, the initial magnetized states of Co 
particles significantly change. The distributions of H, and Hex values therefore are greatly affected, 
which are critical for what the spectra will look like. A simulated spectrum of intrinsic permeability 
in this work is a result of weighted average of many LLG type natural resonances. In addition, by 
controlling the remanence states or particle shape distribution (AN), we can easily control the 


distribution of He, therefore control the permeability spectra [15]. 


0 10 20 30 


40 50 60 70 80 90 100 
f (GHz) 


Figure 6. Spectra of composites only considering Co phase with different volume fractions (a) Co, 


vol.9.63%. (b) Co, vol. 37.51%. 
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When both Co and Fe nanoparticles coexist (see Fig. 2), and are simulated for their microwave 
permeability spectra, which are shown in Fig. 7 and marked as (u’'(AB), u”(AB)). Furthermore, the 
spectra with the mark “(A+B)” indicate that we simply add up the permeability spectra of Fe phase 
and Co phase, as shown Fig.4 and Fig. 6. It is found that when the volume fraction of Co is much 
smaller, the difference between the “(AB)” spectra and the “(A+B)” spectra is not significant. When 
the difference in volume fractions of Fe and Co is not large, the discrepancy of the “(AB)” and 
“(A+B)” spectra is then obvious, as shown in Fig.7b. It is believed this is because the interaction 


among Fe particles and Co particles in Fig. 1b gets stronger. 


SWAB) 


Figure 7. (a) Fe: 90.37% + Co: 9.63%; (b) Fe: 62.49% + Co: 37.51%. 


Such interaction scenario among particles universally exists in electromagnetic composites, which 
can strongly affect the static magnetic parameters, such as coercivity, remanence. In our previous 
work on the permeability spectra of Fe nanowire array, it was showed that when the spacing between 
nanowires is small enough. Localized nanowires showed different coercive fields indicating the 
strong interactions between the nano entities [16]. What is more we can learn from this finding? Let 
us suppose that we have an unchanged particle distribution (exactly same particles) and a fixed 
volume fraction for either Fe or Co particles in their spontaneously magnetized states (i. e. 


unsaturated). At each time, we just mix these same particles differently together and measure the 
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permeability spectra of the composites. Since it is absolutely impossible to ensure the surroundings 
of each particle in these composites unchanged, so do the interactions of magnetic particles. Let alone 
the variable particle distribution of each sampling from a pile of particles powder, it will lead to 
unexpected spectra shapes. Therefore, we can claim that it is impossible to experimentally reproduce 
the permeability spectra. A typical application of electromagnetic materials is to develop 
electromagnetic wave absorbing composites containing randomly oriented particles. For this 
application, the absorbing properties are critically dependent on the permeability spectra. As 
discussed above on the unexpected permeability spectra, it is therefore reasonable to believe that the 
measured performances of electromagnetic wave absorbing composites are uncertain. We prefer to 
name it as an “Uncertainty Principle”. It should be noticed that when we evaluate whether a specific 
magnetic material is good for absorbing electromagnetic wave, we cannot make judgements based 
on several measurements. On the contrary, we need to make a lot of measurements (i.e., statistically 
sampling) to get information about what is the probability to obtain the “good” properties. 
Furthermore, if we also take account of the variation of particle distribution in different measurement 
samples containing the same magnetic fillings. Even though these fillings were prepared under an 
exact same preparation condition (typically, ball milling). The above viewpoint about the uncertainty 
of permeability spectra is more reasonable. Although there are many mixing laws to study the 
permeability spectra of composite, they ignore the fact that the “uncertainty” shown in this paper 
exist from sample to sample, such as the statistical distributions of particle shapes (average, standard 


deviation, median, mode), particle agglomeration and initial magnetization configuration patterns. 


3.2. Effect of damping constants 


For single LLG natural resonance, the damping constant (a) decides how fast the magnetization 
(M) comes back to its equilibrium position during precession, and is positively related to the 
magnitude of high frequency magnetic loss. Actually, when there are plenty of LLG natural 
resonances in a material, the term of “damping constant” is misleading. The intrinsic origins of a 
include the spin-orbit coupling and the scattering by the magnetic impurities. The intrinsic a is NOT 
constant and space-dependent and strongly related to the nonuniform precessions (as shown in this 
work) [17-20]. Besides, the inhomogeneity of microstructure gives rise to the distribution of damping 
constant (known as the extrinsic a) [21-23]. However, there is only one damping constant used in 
LLG equation, which does not agree with the reality. Here, simulations are to show the impact of 
different damping constants on the high frequency permeability spectra, see Fig. 8. Firstly, the 
composite with two phases (Fe: volume fraction = 90.37 % and Co phase: volume fraction = 9.63 %) 
is assumed to have a homogenous damping constant (a=0.5). The spectra of composite are simulated 
and marked with “(AB)”. Secondly, the spectrum of each phase is separately simulated with different 
damping constants (for example, Fe: a=0.3; Co: o=0.7) and added up to get the permeability spectra 
of composite (marked with “A+B”). The “AB” spectra and “A+B” spectra agree badly, which means 


that the inhomogeneous damping constant matters. 
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Figure 8. Spectra of composites if two phases having different damping constants (Fe: a=0.3; Co: 


a=0.7), and two phases having same damping constant (a=0.5, see spectra marked “AB”). 


In addition, with the advances of modern magnetism, especially spintronics, it is found that the 
damping behaviors of LLG natural resonance can be affected by the spin transfer torque (STT). We 
have reported that the imaginary parts of permeability can be negative (u”<0) under the impact of 
STT effect [3]. Besides, the microwave magnetic losses can be suppressed or enhanced by STT effect 
by affecting the effective damping constant. Accordingly, the permeability spectra also can be 


actively controlled by the STT effect [24]. 


3.3. Effect of the remanent states (Memory Effect) 


In practice, the high frequency properties of electromagnetic composites containing magnetic 
fillings are often deteriorated due to the decay of their magnetic properties. The equilibrium 
magnetized states (such as the ones discussed above) are easily affected by the external magnetic 
field, or fluctuant temperature. Accordingly, the permeability spectra will be unstable. Can we 
recover the permeability spectra? Technically, it is feasible as long as we magnetize the materials 
again to establish the original remanent state. Here, we demonstrate the dependence of permeability 
spectra on the remanent states. Firstly, we start to magnetize the composite (Fe: vol. 62.49%; Co: vol. 
37.51%) along the “x” axis with different maximum magnetic fields (see A, B and C point where A 
= 50 mT, B = 80 mT and C = 20 mT) from the initial fully demagnetization state (“O” point in Fig. 
9a); Secondly, the magnetizing field is slowly reduced to zero, and the sample will be in different 
remanent states, see A’, B’points. The composite is the exactly same as the one in Fig. 2b but with 
a slightly different initial magnetization state. The dynamic permeability spectra are simulated as 
before. It is found that same remanent states will have the same permeability spectra. For instant, 
different maximum magnetizing fields (A and C) have same remanent state (A’) which accordingly 
have same dynamic permeability spectra. We call this phenomenon: “Memory Effect”. The word 
“memory” means that if the remanent state is reestablished, the permeability spectra will reappear 
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just like from the “memory” of the magnetic material. The magnetization vectorial distributions of 
these remanent states are obviously different, see Fig. 9(b, c, d). Different remanent states (O, A’ and 
B’) give rise to different permeability spectra, see Fig. 10. Since the difference of remanent state “O” 
and A’ is not sharp, their permeability spectra are similar. But, the spectra of the state B’ and the state 
O are distinct. It can be inferred from these results that any factors strongly influencing the remanent 
state of sample will give rise to different dynamic magnetic response. Hence, when we compare the 
results of an exactly same material from different groups, we should bear in mind that their remanent 
states might be different. Also, it should be noticed that the abscissa axis (frequency) adopts the 
logarithm scale in Figs.10. The differences in permeability spectra between Fig. 10a and Fig.7b (see 
spectra labelled “AB”) are due to their differently initial magnetization configurations although they 


are the exactly same composites. 
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Figure 9. Magnetizing curves and different remanent states 
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Figure 10. dynamic permeability spectra for different remanent states. 


In addition, it is interesting to find that the well-known Snoek’s law does not hold in Fig. 10 a and 
10.c. Snoek’s law tells that the loss peak frequency (f, ) is inverse to the initial permeability (44), 


which means that we can increase the loss peak frequency at the expense of permeability, and vice 
16 / 19 


versa. By controlling the remanent states, the distributions of components of He (i. e., Ha, Hex) can 


be changed [15]. Accordingly, the permeability spectra can be significantly changed. 


4. Conclusions 


For composites containing different volume fractions of Fe and Co particles, their intrinsic high 
frequency permeability spectra have been simulated. A law (called “Mghan’s law”) is proposed to 
explain the spreading out of permeability spectra. It is believed that the local effective magnetic field, 
interaction of magnetic particles, inhomogeneous damping constants and the initial remanent states 
strongly impact on the spectra. The plenty of LLG natural resonances are due to the inhomogeneous 
and local magnetic fields. An uncertainty principle is believed to hold for the electromagnetic 
properties (such as EM wave absorption properties) of composites, which use the magnetic particles 


as the working fillings. Memory effect of permeability is also reported. 
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